Realization of a General Three-Qubit Quantum Gate 



Pasadena, C A 91109-8099 



Farrokh Vatan and Colin P. Williams 
i Jet Propulsion Laboratory 

■ CaUfomia Institute of Technology 
^ ■ 4800 Oak Grove Drive 

a 
o 

^ ; Abstract 

■ We prove that a generic three-qubit quantum logic gate can be implemented using at most 98 one- 
00 . qubit rotations about the y- and z-axes and 40 CNOT gates, beating an earlier bound of 64 CNOT 
^2 ! gates. 

o 

■ 1 Introduction 

O 



Currently, the quantum circuit model is the dominant paradigm for describing the complexity of imple- 
menting a desired quantum computation |1|. By decomposing an n-qubit quantum computation into a 
sequence of one-qubit and two-qubit quantum gate operations, one can characterize the complexity of a 



^ , . quantum computation via the depth of the minimal quantum circuit that implements it. Equivalently, a 

physicist can gain insight into the complexity of performing certain desired manipulations on multi-partite 
k3( ' states. Clearly, any such statement of minimality must always be made with respect to a particular uni- 

■ versal family of quantum logic gates. The most widely used family today is the set of all one-qubit gates 

^ and CNOT gates fljll . However, many equally good universal gate families exist, and the choice of which 

family to use is determined by whichever gates are the easiest to implement in a chosen physical scheme 
for quantum computation. For example, optics-based, superconductor-based, and spin-based quantum 
computers would most likely use CNOT, iSWAP and VSWAP gates, respectively, as their preferred two- 
qubit entangling operation. However, as there are simple relationships between one entangling gate and 
another, insights into circuit complexity based on the {one-qubit, CNOT} family are readily translated 
into the other famihes. Consequently, in this paper we will concentrate on characterizing the complexity 
for achieving arbitrary three-qubit quantum computations using one-qubit gates and CNOT gates. Further- 
more, rather than allowing any one-qubit gate as a primitive, we will only allow the use of phase gates and 
rotations about the y- and z-axes. 

It is known that a maximally general n-qubit quantum logic operation can be implemented using 
0(2^") two-qubit gates l26l . Recently, this estimate has been made more precise for the case of two-qubit 
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quantum computations. Specifically, Vidal and Dawson, Vatan and Williams and Shende, Markov and 
Bullock all proved that a maximally general two-qubit gate can be achieved in a quantum circuit that uses 
at most three CNOT gates. This result remains valid if we replace CNOT with any other maximally entan- 
gling two-qubit gate, such as fSWAP which is the more natural entangling operation in superconductor- 
based quantum computers. Moreover, if the available two-qubit operation, U, is less than maximally 
entangling, than any two-qubit gate can be achieved with at most six calls to U l27l . 

Unfortunately, the aforementioned bounds on the complexity of achieving maximally general two- 
qubit quantum gates have not yet led to similarly tight results for three-qubit gates. Three-qubit states 
are especially interesting because they possess much richer entanglement properties than two-qubit states 
I24II25I I22I. In particular, they ought to allow physicists to gain a much deeper insight into the distinction 
between nonlocality and entanglement fl^ . In order to investigate these states experimentally it will be 
necessary to understand how to best achieve arbitrary three-qubit gate operations. Such an understanding 
might stimulate the development of new quantum information processing protocols that rely upon tri- 
partite entangled states (23). 

To date, there have been relatively few results on the complexity of implementing general three-qubit 
operations. We refer to the most recent paper 1261 . where it is shown, besides the other things, that a 
three-qubit operation can be implemented by using 136 one-qubit gates and 64 CNOT gates. We improve 
this result by providing a computation that utilizes only 98 one-qubit gates and 40 CNOT gates. 

Notation 

We use the following simple notations for one-qubit gates: 

I = ]l2, X = ax, Y = ay, Z = a^,, H = ^ 

For representing the tensor products of several one-qubit operations we apply the following natural nota- 
tion: 

XXX = X®X®X, YYX = y O y O X, etc. 
The following two three-qubit operations play the key rule in our construction: 

N{a,b,c) =exp{i{aXXZ + bYYZ + cZZZ)), (1) 
M{a,b,c,d) =exp{i{aXXX + bYYX + cZZX + dIIX)). (2) 

2 Construction 

The key element of our construction is the decomposition of a three-qubit unitary gate is the decomposition 
of such operations provided by Khaneja and Glaser 1 17 1. This is a general result which recursively reduces 
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the computation of an arbitrary unitary operation on n qubits, U G U(2"), into of a sequence of operations 
on n — 1 qubits and handful "core" operations on n qubits. For example, in the case of n = 2 qubits, for 
U G U(4) we have the decomposition U = Ki AK2, where Ki and K2 belong to U(2) U(2), i.e., the 
space of operations on each qubit separately, and A is a two-qubit operation of the form ^"■^^+I'YY+cZZ ^ 
This decomposition is, in fact, the cornerstone of the recent results on realization of two-qubit operations 
imiUEl. Here we utilize the special form of this decomposition for the case of n = 3 qubits. 

2.1 Khaneja-Glaser decomposition 

In (m it is shown that every unitary operation on three qubits can be decomposed as Figure Q where 
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Figure 1: Khaneja-Glaser decomposition of a three-qubit unitary operation. 

Aj e U(4), Bj £ U(2), Uj = N{aj,bj,Cj) and V = M{a, b, c, d) (see definitions O and (EJi). We utihze 
the optimal construction of |61 for computing the two-qubit operations Aj's. To complete our construction, 
we need computation of the three-qubit operations N{a, b, c) and M{a, 6, c, d). In the next sections we 
provide such computations. 



2.2 Computing N{a,b,c) 

Since the operations XX Z, YYZ, and ZZZ are mutually commuting, we can write 

N{a, b, c) = exp {i{a XX Z + b YYZ)) ■ exp {i c ZZZ) . 
Therefore, we break down the computation of N{a, b, c) to computing unitary operations 
Ni{a,b) = exp {i{aXXZ + bYYZ)) and iVsW = exp [icZZZ). 
First, to compute Ni{a, b), we introduce the following block-diagonal matrix 

P{a,b) = 



Pi 



P2 



(3) 



3 



Pi 



where 

' cos(a — 6) i sin(a — 6) 
I i sin(a — 6) cos(a — 6) 

Then we have the decomposition of FigurelUfor Ni{a, h). 
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Figure 2: Decomposition of the unitary operation Ni{a, b). 

Now to compute the operator P(a, b), we use the technique of Song-Williams fSl for decomposing 
block-diagonal unitary matrices. The result is the circuit of Figure|3l 
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Figure 3: Decomposition of the unitary operation P{a,b). 



We combine the Figures |2 and 15] to obtain a decomposition for Ni{a, b). In this process, we commute 
the SWAP gate with the sequence {l2'^H) ■ CNOT • (12 ® , and next we ehminate both SWAP gates 
by interchanging the rule of the second and third qubits. The result is the circuit of Figure|3 
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Figure 4: Decomposition of the unitary operation Ni{a, b). 



Finally, the circuit of Figure 15] computes the operator N2{c). 

Now by combining the circuits for Ni{a, b) and N2{c) we obtain the circuit of Figure |6lfor computing 
the unitary operation N{a, b, c). 
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Figure 5 : Decomposition of the unitary operation (c) 
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Figure 6: Decomposition of the unitary operation N{a, b, c). 

2.3 Computing M{a, b, c, d) 

Like the previous case, commutativity implies that 

M{a,b,c,d) =exp {i{a XXX + bYYX)) • exp {i{cZZX)) • exp {i{dIIX)). 

So we break down the decomposition of M(a, 6, c, d) to the task of computing the following unitary 
operations: 

Mi(a,6) = exp {i{aXXX + bYYX)), M2(c) = exp {icZZX), and M^{d) = exp {idllX). 

At first step, we have decomposition of Figure for the unitary operation Mi (a, b), where the operation 
P is the same operation defined by 



Figure 7: Decomposition of the unitary operation Mi (a, b). 

Then using the circuit of Figure |3l and eliminating the SWAP gates, we obtain the circuit of Figure [8l 
for computing Mi{a, b). 

The unitary operation M2(c) is simply decomposed as the circuit of Figure |5| 
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Figure 8: Decomposition of the unitary operation Mi (a, b). 
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Figure 9: Decomposition of the unitary operation M2(c). 
Finally, M3((i) = 14 (g) Q, where 



cos((i) i sm{d) 
i sm{d) cos((i) 



i?,(f)-iij,(2d)-i?,(-f). 



Therefore, by putting together all these pieces, we find the circuit of Figure^] for the unitary operation 
M{a, b, c, d) (here we use the identity Q ■ H = H ■ Rz{2d)). 
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Figure 10: Decomposition of the unitary operation M{a,b,c,d), where A = Ry{2a), B = Ry{—2b), 
C = Rz{2c), D = R,{2d), and S = Rz{^). 



Theorem 1 Every unitary operation on three-qubits unitary can be computed by a circuit consisting of at 
most 98 one-qubit phase gates and rotations about the y- and z-axes and 40 CNOT gates. 

Proof. To count the number of gates, we substitute the circuits of Figures |6land[lOlin Fi gure[2 But 
before we start counting the number of gates right away, we should consider the possible "cancellations". 
First note that the first (left) three gates of the circuit of Figure |6l will be "absorbed" by their neighboring 
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gates Ai, A3, Bi, and B^. Also the (only) Rz['^) gate of the same circuit commutes with the gates on its 
right and will be "absorbed" by the gates A2 and A^. Moreover, the last (right) three gates of the circuit 
of Figure will be "absorbed" by their neighboring A3 and ^3 gates; and the left hand side gate S* 
commutes with the three CNOT gates on its left and will be "absorbed" by the gate B2. Therefore, each of 
the operations Ui and U2 contributes 5 one-qubit gates and 9 CNOT gates, and the operation V contributes 
6 one-qubit gate and 10 CNOT gates. 

Finally, we utilize the identity H = az ■ Ry{^), and the fact that every two-qubit unitary operation can 
be computed with at most 15 elementary one-qubit gates and 3 CNOT gates O. I 
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